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Composition of Primary Cosmic Rays at Energies 10 15 -10' 6 eV 
Inferred from Mt.Fuji Emulsion Chamber Experiment 
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The TEy spectrum in 10 2 ~ 5 x 10 3 TeV observed at Mt.Fuji ( 
-v 3725 m a.s.l.) strongly suggests that the flux of primary protons 
in 10 I5 M0 16 eV is lower by a factor of 2-^3 than a simple 
extrapolation from lower energies; the integral proton spectrum 
tends to be steeper than FT 1 7 around 1 0 14 eV and the spectral index 
becomes ^2.0 around 10 15 eV. If the total flux of primary 
particles has no steepening up to MO 15 eV, then the fraction of 
primary protons to the total flux should be ~ 20 % in contrast to 
^45% at lower energies. 


1. Introduction. As well known, the integral energy spectrum of all 
primary particles steepens to FT 2 0 over ^5xlO' 5 eV in contrast to FT 1,6 at 
lower energies. This "knee” of the spectrum may come from 1) a change of 
the acceleration mechanism resulting in a steeper source spectrum and/or 
2) a breakdown of magnetic trapping in the galaxy. Study on the 
composition around the "knee” is an important key to resolving the 
problems of the origin and propagation of cosmic rays. However, the 
definite result has not yet been obtained because of the too low primary 
flux for the direct observation and too poor response of the air shower 
observation to the primary mass. 

The emulsion chambers ( EC ) at Mt.Fuji have been extended to the 
scale of 1000 m 2 x year, and a few hundred events produced by primaries 
with energy 1 0 l5 -^- 1 0 1 6 eV have been accumulated 1 ^ . Meanwhile, accelerators 
has clarified the feature of nuclear interactions up to 150 TeV (in lab. 
energy). These circumstances enable us to investigate the primary 
composition around the "knee” by analysing the EC events through 
simulations based on a reliable model of nuclear interactions. 


2. y-ray Family. A group of high energy particles (mostly y-rays with 
energy > a few TeV) produced by successive collisions of a primary 
particle with air nuclei is called a family. Particles which penetrate 
deep in the atmosphere can produce families effectively, otherwise 
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produced particles are usually absorbed in the atmosphere. Therefore 
protons are most efficient in producing families; the average energy of 
protons responsible for y-ray families with EEy>100 TeV is -^2xl0 13 eV 
while that of irons is almost one order of magnitude higher. 

It has been found that the fraction of proton-induced families in 
all families in E£ y >100 TeV is ~ 70% even if protons constitute only ~ 
10% of all particles around 10 15 eV. Therefore, if we convert the family 
intensity into the primary flux assuming that all families are produced 
by protons, we obtain an upper bound of the proton flux within an 
accuracy of 30%. 


3, Monte-Carlo Simulation. The Monte-Carlo method has been used to 
simulate the family phenomena and to deduce the primary intensity by 
comparing the obtained family flux with the experimental one. Important 
factors in the simulation are the cross-section of nuclear interactions, 
the leading particle spectrum and the secondary particles spectrum 
(especially, the degree of scale breaking in the fragmentation region). 
These are summarized in the following. 


Cross_ Section Using the 
data of pp and pp interactions up 
to 150 TeV, the mean free path of 
the N-Air collision is expressed as 
A=Aox£T 5 , where E is in TeV, 6 is 
0.04 ~ 0.06 and Ao=82 g/cm 2 . 

Elaslisity Two cases of 
the elasticity distribution are 
assumed as shown in Fig.l. Case A 
disregards the difference between 
protons and neutrons and gives the 
average elasticity of 0.44. In case 
B is considered the absence of the 
diffraction peak in charge exchange 
processes, for which we refer to 
p-C and p-Be data 3) at 100 - 400 
GeV. The average elasticity in a 
non-charge exchange process like 
p+Air—p+X is 0.45 while that of 
p+Air— n+X is 0.28. With a charge 
exchange probability of 28 %, the 
average elasticity in case B 
becomes 0.40. 



ELASTICITY ( X ) 

Fig.l Elasticity distribution. 


E nergy. distribution _o£ secon d ary gartiolss The QCD-jet 2) 

automatically introduces a mild scale braking in the fragmentation as 
shown in Fig. 2a. The rapidity distribution shown in Fig. 2b is compatible 
with the CERN-SPS data. To see the effect of a strong violation of 
scaling in the fragmentation region we employed a fireball model of which 
the multiplicity grows as E 1 / 4 . 

In Table 1 , the models are characterized by the parameters discussed 
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above. The calculated flux of y-ray families at E£y=100 TeV for each 
model is listed by normalizing the flux to model QA6. 




Fig. 2a x-distribution of Fig. 2b y-distribution 

secondary pions in p-Air corresponding to Fig.l. Black 

collisions at 10 4 * * * * * * * 12 * eV and 1 O 16 dots show data by the CERN-SPS. 

eV. A mild-violation of scaling Correction due to Air-target 
is seen in the model with QCD effect is included, 
jet. 


Table 1. Interaction model and the relative flux of y- ray family. 


model 

x distribution 

elasticity 

m.f .p. (5) 

flux ratio 

QA4 

QCD 

case A 

0.04 

1.32 

QA6 

QCD 

case A 

0.06 

1 

QB6 

QCD 

case B 

0.06 

0.72 

FA4 

Fireball 

case A 

0.04 

0.41 


4. Proton Spectrum in 10 15 -^10 16 eV. Under the assumption that all 

y-ray families are produced by protons, an expected spectrum of protons 

is deduced as shown in Fig. 3. The upper bound is given by model QB6 and 

the lower one by model QA4. As mentioned already, a maximum of 30 % 

decrease might take place in the expected flux (due to the heavy primary 

contribution). Even if we take the experimental errors (both statistical 

and systematic) into account, we are forced to conclude that the proton 

spectrum should steepen around 10 t4 eV, and the ratio of protons to all 

particles in 10’ 5 ~ 10' 6 eV should be lower by a factor of 2~ 3 than that 

at lower energies. 
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E/Nucleus ( TeV ) 

Fig. 3 Expected primary proton spectrum. 

5, Discussion. By an analysis of air shower data, some authors 
reported a proton-dominant composition^ . In their estimation, a scaling 
violation in the central region are included to be consistent with the 
CERN-SPS data while the x-distribution in the fragmentation region is not 
properly considered, which is most effective to the cosmic-ray 
propagation in the atmosphere. Even if we assume a strong scaling 
violation (FA4) which is unrealistic already at the accelerator energies 
as shown in Fig. 2b, the expected proton flux never exceeds the simple 
extrapolation from low energies. 

The general (i.e., all observed particles independent of family 
correlation) hadron and y spectra observed at Mt.Fuji also supports the 
present picture 5 . In Ref .6, we estimated the composition around the 
’’knee" by the use of the “rigidity cut-off model”. The proton spectrum 
there is consistent with the present estimation. 
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